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Abstract 

Using a new approach, we propose an analog of the Fizeau effect for massive and 
massless particles in an effective optical medium derived from the static, spherically 
symmetric gravitational field. The medium is naturally perceived as a dispersive 
medium by matter de Broglie waves. Several Fresnel drag coefficients are worked 
out, with appropriate interpretations of the wavelengths used. In two cases, it turns 
out that the coefficients become independent of the wavelength even if the equiva- 
lent medium itself is dispersive. A few conceptual issues are also addressed in the 
process of derivation. It is shown that some of our results complement recent works 
dealing with real fluid or optical black holes. 

PACS number (s): 03.65.Bz, 42.15.- i, 04.20.-q 



1 Introduction and reappraisals 

The historical Fizeau effect for light in moving media has been reconsidered by several 
authors [1-5] in recent times. We shall consider it here in the context of static, isotropic 
gravity. To our knowledge, such an investigation has not been undertaken before. We 
deemed it worthwhile to examine how an old effect would look like in a new theoretical 
model and what conceptual issues are involved. However, we must make clear at the outset 
that the only quantity to be borrowed from general relativity is the effective refractive 
index. The rest of the analysis is special relativistic (see Sec. 3). 
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In the literature, generally, the Fizeau effect is considered in connection with its close 
relative, the analog of the Aharonov-Bohm (AB) effect in real material medium. Sev- 
eral interesting results have followed from these analyses. For instance, Leonhardt and 
Piwnicki [6] have shown that a nonuniformly moving medium appears to light as an ef- 
fective gravitational field for which the curvature scalar is nonzero. They also show how 
light propagation at large distances around a vortex core shows Aharonov-Bohm (AB) 
effect and at shorter distances resemble propagation around what are termed as optical 
black holes. Berry et al [7] demonstrated the AB effect with water waves and Roux et 
al [8] observed it for acoustical waves in classical media. The curved space analogy has 
been predicted for fluids and superfluids as well [9]. The spirit of the present work, in 
some sense, is in a direction that is reverse to the idea of the above curved space analogy. 
That is, our interest is to calculate the Fizeau type effect for both massless and massive 
particles in a static, spherically symmetric gravity field hut portraying it as an effective 
optical medium. In the process, we shall also see the extent to which the curved space 
analogy compare with the results derived in a genuine gravitational field. For simplicity, 
we shall assume only uniform motion of our effective medium resulting from the relative 
motion between the gravitating source and the observer. 

An outline of the Fizeau effect is this: Consider a tube through which a fiuid with 
a refractive index n is fiowing with velocity V. Then, let light pass through the tube 
parallel to its axis. In the comoving frame of the water, the speed of light is v'{= Co/n), 
but in the frame in which the water appears to be flowing, the speed of light has been 
found to be: 

. = „'±(i-Jv + o(v^)«^±(i--L)v. (1) 

where cq is the speed of light in vacuum. The quantity (1 — is called the Fresnel 
drag as he was the first to predict it theoretically. Obviously, the resultant speed v does 
not conform to the Galilean law of addition of velocities v' ±V. The effect, after it was 
experimentally observed by Fizeau in 1851, was regarded as an empirical fact awaiting 
a correct theoretical interpretation. It came only after the advent of Einstein's special 
theory of relativity in 1905. It has since been realized that the Fizeau effect symbolizes 
only a first order approximation of the exact one dimensional special relativistic velocity 
addition law (VAL) derived from Lorentz spacetime transformations. Originally, Fizeau 
did not consider dispersion but nowadays it is recognized that the effect also contains a 
term due to the effect of dispersion. 

In our investigation, we shall adopt an approach involving quantum mechanics, gen- 
eral and special relativity using the method of what is known as the optical-mechanical 
analogy. The historical and fundamental role of the analogy in the development of modern 
theoretical physics need not be emphasized. Apart from the crucial role it played in the 
development of quantum mechanics, especially in the de Broglie wave-particle duality, it 
provides an excellent tool that enables one to visualize the problems of geometrical optics 
as problems of classical mechanics and vice versa. 

In a series of papers [10], it has been shown that the optical-mechanical analogy can 
be recast into a familiar form that allows one to envisage the mechanical particle equation 
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as a geometrical optical ray equation and the latter as a Newtonian "F = ma" equation: 



§-v(in^cS). (2) 

dA = (3) 

where r = {x, y, z) or (r, y?) , V is the gradient operator, n is the index of refraction, not 
necessarily constant, and A was originally called the stepping parameter but it could also 
be identified as the optical action and related to several other physical parameters. Many 
illustrations in ordinary gradient index optics demonstrated the validity of the Eqs.(2) 
and (3) and their usefulness heuristic tool. 

An interesting turn in the direction of investigation is signaled by the introduction 
of general relativity [11]. Exact equations for light propagation in the static, spherically 
symmetric field of Schwarzschild gravity do indeed follow from the Eqs. (2) and (3) when 
an appropriate gravitational index of refraction n{r) is employed. The analysis also brings 
forth the distinct but complementary roles played by the optical action A and coordinate 
time t. To see this, note that the first integral of Eq.(2) is 



dr 



dA 



or equivalently, using Eq.(3), 

dr 



nco, (4) 



dt 

However, the force laws have changed thereby. In Eq.(4), the "potential" is |n^CQ while in 



Eq.(5), the potential is On eliminating A from Eq.(4) or t from Eq.(5), we would 

therefore obtain two path equations for light on a plane, but only the former, not the 
latter, gives the right answer. On the other hand, Eq.(5) gives the correct equation for 
the Shapiro time delay At, while AA from Eq.(4) does not. A deeper understanding of 
the parameter A is still awaited. 

Our basic strategy is to regard the gravity field as an effective refractive optical medium 
imposed on a fictitious Minkowski space so that Lorentz transformations can be used to 
relate two relatively moving observers in that space. (Note that we are not talking of a 
division of the metric tensor into two parts, but rather of a scalar field placed upon a flat 
space. For more discussion, see Sec. 3). This is just an intermediate exercise. The final 
outcome has to be translated back into the actually observable quantities in a gravity 
field. The idea that a gravity field could be formally equivalent to a refractive medium 
with respect to optical propagation is not new. It goes all the way back to Eddington 
[12] who was the first to advance the expression of a gravitational refractive index in 
an approximate form. It was used later, in varying degrees, by several other researchers 
[13,14] in the investigation of specific problems. But none of the works really focused on 
how the exact general relativistic equations of trajectories, frequency shifts or Shapiro time 
delay for massless particles could be obtained in that equivalent medium. The motion 
of massive particles was not addressed at all. The extension of the work in Ref.[ll] that 
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includes also the massive particle motion now exists [15,16]: A suitably modified index of 
refraction together with the "F = ma" formulation immediately reproduce all the desired 
exact equations in the static, isotropic gravity field. The method has been applied very 
successfully to Friedmann cosmologies as well that yielded some new interesting insights. 
All the above systematic developments amply indicate the usefulness of the concept of an 
effective gravitational index of refraction. By way of a further extension, the index has 
been calculated also for a more general class of rotating metrics [17] . A new and significant 
development has come in the shape of a most recent formulation [18] of a single set of 
unified optical-mechanical equations that allow easy introduction of quantum relations 
into it. As a consequence, one then finds that massive de Broglie waves necessarily 
perceive the gravity field as a dispersive optical medium. 

In this paper, our basic aim hinges around calculating the consequences arising out 
of this dispersion in the form of what may be termed as the gravitational Frcsncl drag, 
dispersion included. There are several spin offs: It will be demonstrated that, in the 
comoving frame, the expressions for the Lagrangian and the dispersion relation are similar 
to those obtained by Leonhardt and Piwnicki [6] in the context of real media. These 
similarities provide a direct extension of these expressions in a realistic gravity field. It 
will also be evident that the conditions for optical black holes [6,19,20] are naturally met 
in the equivalent medium, irrespective of whether one considers light or massive de Broglie 
waves. 

The paper is organized as follows: Sec. 2 contains a brief survey of the basic equations 
that will be used throughout the paper. Conceptual justifications for the adopted pro- 
cedure appear in Sec. 3 . In Sections 4-6, the gravitational Fresnel drags are calculated 
for different choices of the wavelengths. Sec. 7 contains a brief discussion of operational 
definitions. In Sec. 8, we demonstrate how the results dealing with a real fluid medium 
compare with those in a genuine gravity field considered in this paper. Finally, in Sec. 9, 
we summarize and add some remarks. 



Consider a static, spherically symmetric, but not necessarily vacuum, solution of general 
relativity written in isotropic coordinates 



where f2 and $ are the solutions of Einstein's field equations. Many metrics of phys- 
ical interest can be put into this isotropic form including the experimentally verified 
Schwarzschild metric. The coordinate speed of light c(r*) is determined by the condition 
that the geodesic be null (ds^ = ): 



We take leave from the metric approach at this point and define the effective index of 
refraction for light in the gravitational field as 



2 Basic equations 




(6) 




(7) 



-1 



(8) 
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We shall omit further details here that can be found in Ref.[18], but only state the 
results to be used in this paper. The first step in the direction of introducing quantum 
mechanics in a semiclassical way is to have a single refractive index N and a single set of 
equations that should be valid for both massless and massive particles. The result is: 



dr' 



V -N^Cq ' (light and particles) (9) 



dA 



Nco, (light and particles) (10) 



where, once again, it is the same A, satisfying dA ~ that appears even for massive 
particle trajectories. It looks as if the action has a foot in the wave regime and a foot in the 
particle regime. The second step involves the introduction of the Planck relation H' = huj' 
and the de Broglie relation p' = hk' = h/X', where h = 27rh, in the expression for N. 
As usual, H' and p' are the total energy and momentum respectively, uj'{= 2ttv') and A' 
are the coordinate frequency and the wavelength of the de Broglie waves. The physically 
measurable corresponding proper quantities are a)' = uj' /Q. and A' = A'/$ respectively. 
The third step finally gives the desired index of refraction N of the dispersive medium 
due to massive de Broglie waves: 



N{f^u')^n{n\ll- '^"f^P ^ (11) 



where m is the rest mass of the test particle. One may also rewrite N as 

H' Co ' ^ ' 

where v' is the (unobservable) coordinate speed of the classical particle in the medium. 
Also it follows that 

-' = 2..', A'=^. (13) 
Using Eq.(13), N can be rewritten in a more transparent form: 

iV=^£) , (14) 



where Ac = h/mco is the Compton wavelength of the particle. Clearly, for light, m — 
0,N — n, and one recovers Eqs.(2) and (3) from Eqs.(9) and (10) respectively. That 
is, light waves do not perceive the effective medium as dispersive. However, for m ^ 0, 
dispersion seems inevitable, as evidenced from Eqs.(ll) or (14), if quantum relations are 
introduced. 

We shall require also the following: The mass shell constraint is given by [18]: 



n uj — m Cr,^i H — . (15) 
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The phase velocity is 



p' /e' iv " n 

giving the group velocity 



duj' CqN I 
dk' V? 

It should be mentioned that the validity of the expression (11) is established also by the 
WKB analysis of the massive generally covariant Klein-Gordon equation [18]. Moreover, 
the mass shell constraint (15) yields the exact Stodolsky phase [21] in the case of spin-1/2 
Dirac equation in curved spacetime [22]. This last result is extremely interesting. 

With the Eqs.(6)-(17) at hand, we are able to calculate the Fresnel drag factors under 
different scenarios, but, before this, we need to clear up a few relevant concepts. Note 
that all the expressions in this section refer to the comoving frame, that is, the frame 
fixed to the gravitating source. Henceforth, in order to have conformity with notations 
in the literature, all expressions in the comoving frame will be designated by primes and 
those in the relatively moving lab frame will be denoted by unprimed ones. 



3 Conceptual issues 

The following discussion is aimed at providing appropriate interpretations of the quan- 
tities that appear in the various formulations of the Fizeau effect. There are two basic 
ingredients. The first is the VAL. In many works dealing with the effect, the one di- 
mensional VAL, which is valid for point particles, is employed, implicitly or explicitly, 
also for waves propagating with the phase speed co/n. The procedure is to use the one 
dimensional Lorentz transformation equations in the form 

u' = ^{uj-kV), k' = ^{k-Vujcf), (18) 
uj^^{J + k'V), k^^{k' + VJcf), 7 = (1 - 1/Vc2)-V2 (19) 
and obtain a VAL as 

v'p = (vp -V){1- v,V/cl), Up = (.;; + V)il + v'^V/cl), (20) 

where v'^ = uj'/k' = Cq/ti is the phase speed of light in the (primed) comoving frame of 
the medium and Vp = uj/k is the phase speed in the (unprimed) lab frame in which the 
medium appears to be moving with uniform relative velocity V. The phase speed however 
could well exceed cq in many physical configurations where n < 1. 

On the other hand, an a priori prescription that n be greater than unity (making 
Cq/u < Cq) somewhat diminishes the generality of the theory. However, this deficiency 
may not pose any realistic problem in a nondispersive medium. When dispersion is 
involved, the most appropriate quantity to use in the VAL is the group speed duj/dk (that 
involves the knowledge of dn/duj), which simply equals the velocity of the classical point 
particle, rather than the phase speed. As stated before, the original Fizeau experiment 
did not consider any dispersion, the index n was taken to be a true constant, so that the 
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group and the phase velocities coincided precisely to co/n. In general, they are different 
for massive de Broglie waves, as our later equations will reveal. In our calculation of 
the Fizeau effect, the mass shell constraint, Eq.(15), or, by another name, the dispersion 
relation plays a key role: It provides well defined expressions for the group and phase 
velocities. Such types of natural constraints arc unavailable in just any arbitrary medium 
consisting of solids or liquids. In these cases, dispersion is normally introduced by hand. 

An important point should be noted here. In describing the Fizeau experiment with 
ordinary medium (such as water), one takes the background spacetimc to be flat. Such 
Minkowski networks, composed of rods and clocks, are actually unobservable in a gravity 
field due to the universality of gravitational interaction, or, saying more technically, due 
to the principle of equivalence. There does not exist a unique division of the metric 
tensor into a background and a field part. We consider here a different kind of separation 
according to which the gravity field is looked upon as analogous to an optical medium 
imposed upon a flat background spacetime, the index summarizing the nonlinearities 
of the gravity field, as it were. The important point is that the analogy, though intended 
to be only of formal nature, may lead to results that could be testable by experiment (see 
Sec. 9 for a discussion). With this understanding, let us conceive of observers equipped 
with fictitious Minkowski networks and apply, as an intermediate step, the full machinery 
of special relativity in what follows. 

Thus, we take Eq.(19) in the form 

Au = 7( Au;' + Ak'V) , Ak = 7( AA;' + VAu'c^'^) , (21) 
which give the VAL, denoting Vg = as: 

V - < + ^ (22) 

The second ingredient is the special relativistic Doppler shift in one dimension giving 
frequency (or wavelength) transformation between two frames in relative motion. Thus, 
one takes Eq.(19) in the form 

uj^-iUj'{l + k'V/uj'), (23) 

and specifies u' /k'. At this point, let us note that. Cook, Fearn and Milonni [2] have 
considered two possibilities in the context of Fizeau experiment with real media having 
refractive indices n. 

Case 1 

Take u' /k' = cq in Eq.(23). This case has been considered by Synge [3]. That is, take 
the usual Doppler shift formula, which, written in terms of the wavelength, is 



1-^ 

Co 



\ 1 + ^' 

l| CO 



(24) 



The corresponding physical configuration consists of a block of material moving with 
velocity V in an otherwise empty lab frame. The wavelength A' of a light pulse measured 
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by an observer stationed at the interface between the block and the empty space will 
appear to the lab observer as A according to Eq.(24). Inside the block, however, A' is 
assumed to be a constant. The resulting Fresnel drag has been experimentally confirmed 
to a very good accuracy by Sanders and Ezekiel [4]. 

Unfortunately, it is difficult to conceive of a parallel configuration in our problem. The 
entire optical medium can neither be simply put inside a box with a certain boundary 
nor need the wavelength A' be constant throughout the medium. Instead, it is easier to 
consider two relatively moving observers associated with the background empty frame 
who may use Eq.(24). We have to calculate how one observer translates the observations 
of another at a certain point when they happen to pass each other. This is done in Sec. 
4. 

Case 2 

Take uo' /k' = Co/n in the Eq.(19) for k. Cook et al [2] provide the corresponding 
physical configuration in this case. According to Lerche [1], the lab observer can exercise 
two options. Either he/she uses (i) a wavelength A given by the Doppler formula (23) 
but with uj' /k' — cq/u or uses (ii) a vacuum wavelength Aq = 2t:cq/uj. The forms for the 
drag coefficients will be different in the two cases. The parallel options in our case are 
the same, except that we have to use instead of n, so that the Doppler formula reads 

We shall work out both the options in Sec. 5. This particular formula appears to be more 
consistent with our formulation per se as we will be using our own definition of cu'/ k' given 
in Eq.(27). We can also add a third possibility worked out in Sec. 6. This is a special 
feature of the gravitational case we are considering. 

Case 3 

Consider a stationary observer A at a point in the gravity field measuring proper (or 
physical) wavelength A'. He/she also measures the proper velocity of light in his/her 
neighborhood to be just cq. A freely falling observer B at that point, having an instanta- 
neous velocity V relative to A, would measure A according to the options, which, to first 
order, are 

A' A(l + F/co), X' ^X{l + VN/co) Ao = 27rco/a;. (26) 

Note that there is a difference between the present stationary observer and a stationary 
observer associated with the background flat space of Case 1: The group velocities of 
the matter de Broglie waves measured by them are not the same (see below). We now 
proceed to calculate the Fresnel drags successively in all the three cases using the same 
VAL, Eq.(22), but different Doppler formulas, Eqs.(24)-(26). 
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4 Fresnel drag: Case 1 



Suppose that an observer A, equipped with a Minkowski network, is at rest at r = 
in a spherically symmetric medium. He/she measures the coordinate phase and group 
velocities respectively of a massive de Broglie wave packet at r = Tq as, using Eqs.(15): 

k' '""P' N{r'o,X'y ^ ' 

For a light pulse, Vp — v'g = ^ and these are independent of the wavelength A' or wave 
number k'. The same holds for Vg as well. This implies that the trajectories of light rays 
do not depend on the wave properties of light. However, in general, v'py^v'g, as is evident 
from Eqs.(27) and (28). 

Consider another observer B moving in the same radial direction with uniform velocity 
V with respect to A. Then, in the frame of B, identified as the lab observer, the entire 
medium moves uniformly, that is, A becomes the comoving observer. How will B translate 
the observations of A, when their origins coincide at r = 0? To find it out, note that the 
coordinate length Tq will appear to B as 

(29) 

Also, the velocity v'g observed by A will appear to B as Vg given by the special relativistic 
VAL, Eq.(22). We may exphcitly express v'g in terms of (r^,. A') as 

/ Co Co , . 

= TFTZr^ = / . .o =• (30) 




When this expression for v'g is plugged into the right hand side of Eq.(22), one finds the 
answer to the question above: Vg{r'Q, A') is the exact radial group velocity of the de Broglie 
waves to be observed by B. But B uses the Doppler shifted wavelength A instead of A'. 
Then, to first order in (V/cq) ), we get from Eq.(24): 

A' ^ A(l + V/co) = A + AA, r'o^ ro. (31) 

Considering the right hand side of Eq.(30) and writing the denominator as -/V(ro = 
ro. A') = N{X + AA), we get from the Taylor expansion 

dN - XV dN 

From Eqs. (22) and (28), we get, using the above, a redefined index N' such that 
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In other words, in the approximation considered, A^^(A') ^ N'^{\) and we have. 



where 



N{X + AA) 
Fi=(l- 



+ 1 



1 



iV2(A)_ 



V 



A 



X 



7V(A) 

dN 



+ FiV, 



(33) 



(34) 



Ar2(A) ; iV2(A) aA 

is the Fresnel drag we have been looking for. It can be easily verified that the same 
Fi follows also from the ordinary expansion of Vg {tq, A') in Eq.(22) in conjunction with 
Eqs.(21) and (30) under the small velocity approximations, Eq.(31), but the steps as given 
above are the simplest. For hght waves, TV — > n, and one has 



Fi = 1 - 



A dn 

X 



n2(A) ; n2(A) ^A 



(35) 



Interestingly, although the dependence of n on A is not known, the dispersion nonetheless 
follows here as an inheritance from Eq.(34). This is the formula proposed by Synge [3] 
and also experimentally tested [4] with n as the refractive index of the block. 
Using Eqs.(14) and (28), we can have the explicit expression from Eq.(34) as: 



ay 



1 + ^ 



n(ro) 



3/2 • 



(36) 



Note that, in the asymptotic region r — >■ oo, or in the absence of gravity, one has 
n{r) ^1, A — ^ A', so that, from Eq.(28), the group and phase velocities of de Broglie 
waves, as measured by A, respectively are 



Co 



1/2 



< Co, 



Co 



Ac, 



1/2 



> Co 



(37) 



and thus one finds that matter de Broglie waves perceive even the fiat space as a dispersive 
medium with an index of refraction 



Nflat = 



1 + 



A. 



1/2 



(38) 



One recognizes that it is this v'g in Eq.(37), together with 
transformation law: 



Ac 



provides the energy 



H 



mcn 



(39) 
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Then, one recovers the special relativistic mass shell condition. It follows that, in this 
case, the drag measured by B in terms of his/her wavelength A, is 



flat 



X 



1 - 



{0 



1/2 



(40) 



As one can see, Eqs.(37)-(40) arc restatements of the well known special relativistic ex- 
pressions, but only interpreted in a different way. 



5 Fresnel Drag: Case 2 

According to first option (i), the Doppler shift is given by Eq.(25). Thus, we have, to first 
order in (V/cq): 

X' ^ X{l + V/Nco)^X + AX, AA=^. (41) 

N cq 

Then, writing again: N{rQ = ro. A') = N{X + AA), we get from the Taylor expansion 

iV(A + AA) . ^(A) + AAf ^ ^(A)(l + -^f ). (42) 



The resultant group velocity as observed by B, who uses A of Eq.(41), is 

= = ^(l - ^J—\ V 

N'{X') N{X') \ m{X') J 



+^-TfTrT^]v-^ + F2V, (43) 



A^(A + AA) V ^^WJ 
where 

is the drag factor. For light waves, N — n, N = n so that 

„ A 1 \ A dn 

X — . (45) 



This formula was first given by McCrea [23]. Writing exphcitly, we find from Eq.(44), 

F2 = l- (46) 

This coefficient comes out to be independent of A! According to the second option (ii), B 
uses a vacuum wavelength. In this case, the calculations would proceed slightly differently. 
Consider Eq.(18) for lo' instead of Eq.(41). Then, we have, to first order, 

uj' ^uj{l-VN{uj)/co) =uj + Auj, Auj = - '^^^^'^\ (47) 

Co 
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Then, proceeding as before, 



Vg{u) 



Co 



N{uj) 



+ 



_ 1 \ ^ 9N 



V. 



Now B uses the vacuum wavelength to be Aq = 2t:cq/u, so that Eq.(48) gives 



Vg{\o) = 



Co 



A^(Ao) 



where 



For hght waves, we get 



1 



1 \ _ AoAr(Ao) dN 

Ao 



1 



dn 



n2(Ao)/ n(Ao) <9Ao' 



(48) 



(49) 



(50) 



(51) 



This is the expression given by Lerche [1] and Cook et al [2] for Fizeau experiment with 
water with index n. Writing explicitly, we find from Eq.(50), 



1 - 



( 



n2(ro) 



n 2 



(52) 



where Aq = Ao$~^ Thus, so far, corresponding to and iV, we have three Fresnel coeffi- 
cients Fi , F2 and F3 depending on the VAL and the various Doppler shifted wavelengths 
used by B, as considered in the literature. 



6 Fresnel drag: Case 3 

Consider an observer A at rest with respect to the gravitating source at a coordinate 
radial distance r = Tq. He/she will measure proper quantities. The mass shell condition 
would be given by 

h^u'^^m'^cl + clh^k'^, k'^^k', (53) 

so that A measures, in his neighborhood, the proper phase and group velocities of the de 
Broglie waves which are connected by 



v'v' = — ^ = Co, (54) 
" ' k' dk' ° ^ ' 



where, using Eq.(27), 



, uj' ( ~ . / r/N / dw' Co ~ n , ^ 
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Note that, these phase and group velocities are not the same as those measured by A, 
viz., Eqs.(27) and (28), which highhght the difference between the two observers. The 
observer A measures the velocity of light as Vp = v'g — cq since N — n. Consider another 
observer B falling freely in the same radial direction attaining an instantaneous speed V 
at r = Tq. Since the frame in which B is at rest is locally inertial in virtue of the principle 
of equivalence, the speed of light measured by him will also be cq and hence A and B 
can be connected by a Lorentz transformation. Then Vg would appear to S at r = ro as 
■Uggiven by the VAL: 



V'g + V 
1 + ^ 



(56) 



Employing arguments similar to those in Cases 1 and 2, we can straightaway write down 
the corresponding drag coefficients: 

(a) A measures A' and B uses A connected by X' ^ X{1 + V /cq): 



iV2(A) 



JV2(A) 



X 



m 



1+ 



ii)' 



X 



1+ 



T72 



(57) 



(b) A measures A' and B uses A connected by A' A(l + VN/cq): 



1 



^x — = 
N^{X)J iV(A) ^ dX 



(58) 



(c) A measures Cu and B uses Aq connected by Aq = 27rco/a): 



iV2(Ao) 



Ao 



X 



dN 



= 1 



1+ 



X 



Ar3(Ao) " dXo 



1 + 



1+ 



(69) 



where Aq = Ao$ ^. We also see that the radial proper velocity of the classical point 
particle as measured hy A at r = is given by 



v' 

prop 



dr' 



n 



dr' 
It' 



nv 



coord ■ 



(60) 



we find that f^^^p = v'^. For 



Using the definitions: dV = $ ^dr', dr' = Qdt', v'^^^^ 
light, of course, v'^^^ = ^, and v'p^op = % = Cq. The last result is also consistent with the 
fact that ds"^ = c^dr'"^ — dl'"^ = gives Jp- = Cq. For light waves, we find N = 1, so that 
Fi = F2 = = 0. These indicate only the special relativistic invariance of light speed, 
no matter what wavelength B uses. For de Broglie waves, the difference among the drag 
coefficients is evident from Eqs.(57)-(59). 
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7 Operational definitions 



In order to operationally realize the value of Fi in a gravitational field, consider a simple 
thought experiment. Let there be a source in free space that produces de Broglie waves 
with wavelength A'. Then A is known via Eq.(31) which is the wavelength measured by B. 
Let A take this source to any point inside the refractive medium. Then, he will measure 
the same A' as A' = A'^"^ and B will find A = A$~^. The only other quantity is the 
coordinate distance Tq appearing in the refractive index n(ro) and $(ro). The expressions 
for the index is supplied by the metric functions. For instance, in the Reissner- Nordstrom 
field, with G = Co = 1, we have 



1 - 



(M2 - g2 



(61) 



(62) 



where M and Q are the mass and the electric charge. For the Schwarzschild field, we have 
Q = 0, so that 



n{r) 



i + f) 
(i-S)' 



(63) 



If the relative velocity V between A and B is small, V « Cq, we can take Tq 



from Eq.(29). If we consider that both the observers are in a weak gravity field, we can 
take To ~ Tq I, where I is the physically measurable distance from the center of the 
gravitating source to the field point. Then 



n(ro) ~ n{l) ^1 + 



2M 



(64) 



With these inputs, Eq.(36) provides the theoretically predicted value of F-i after the known 
value of the Compton wavelength is plugged in. 

Interesting results are obtained in the case of F2 and F2- One finds that F2 does not 
involve the wavelength at all. This means that a Fizeau type experiment either with 
light or with de Broglie waves would yield the same drag factor, if Eq.(23) is followed 
in conjunction with Eq.(27). In this case, it appears that the wavelength dependence 
introduced by the group velocity is undone by Doppler shift. A similar thing occurs also 
in the case of F2 which is identically zero. 



8 Comparison with real medium 

Starting from the wave equation in a nonuniformly moving fluid with refractive index n, 
Leonhardt and Piwnicki [6] derive the Lagrangian and the Hamiltonian for a light ray 
as observed by a lab observer. From the action principle, they arrive at a completely 
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geometrical picture of ray optics in a moving medium. Light rays arc geodesic lines with 
respect to Gordon's metric, which, in the comoving frame reads 



(65) 



The Lagrangian, Eq.(49) of Ref. [6], they derived for a light particle in the lab frame, is 

(66) 



L — —mco^ 



-1 tMco- 



u.v 



Co 



where u is the fluid velocity in the lab frame, v (= 2 ) is the velocity of the light 

particle conceived of having a fictitious mass m and 7^ = — ^) .In the comoving 



frame of the fluid element, u — Q so that 



L = -mco X - X 



1 
n 



Co 



(67) 



Consider the Lagrangian for a massive particle in the comoving frame, derived in our 
Ref. [18], viz.. 



L — —mc^Q. 



1/2 



(68) 



where v' is the classical particle coordinate speed. Now note that the metric (65) with n 
as the real medium index, can be obtained formally from Eq.(6) above simply by putting 
^ — 1 and Q, — 1/n. Clearly, the n in Eq.(68) has a different origin: it derives from 
general relativity. Using this value of Q in Eq.(68), one flnds that it is exactly the same 
as Eq.(67). 

The dispersion relation for light in the comoving frame {u = 0) following from Eq.(33) 
of Ref. [6] is 

(69) 
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0. 



This is precisely the same as that following from Eq.(15) with m = for light. 

Interestingly, taking a cue from Eq.(66), we may proceed to write down the Lagrangian 
of the classical particle in the lab frame as follows: 

Our metric, Eq.(6) in the comoving frame can be written down as 



ds^ 



$-2 X 



cldt'^ - df^ 



- 1 ] ddt'^ 



(70) 



To go to the lab frame, we effect a Lorentz transformation. Note that there is a Lorentz 
invariant term in the parenthesis and hence only the last term needs to be transformed. 
Thus, in the lab frame the metric is 



ds'^ = g^ydx^dx"" = $" 



X 



1 



- 1 KK 



(71) 
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where ?7^j. = [cq, —1, —1, —1] , = 7 (^1, — , 7 = (^1 — , and V is the veloc- 

ity of our medium in the lab frame. In the comoving frame, = (1, 0, 0, 0). The action 
is given by 

S = -mco J p,.^^dt = J Ldt 
Defining v^" ^ ^ ^ {l,v), we can find the Lagrangian for a particle in the lab frame, 



L = -mco$ ^ X 

The dispersion relation (or, the Hamiltonian) in the lab frame can also be obtained by a 
Lorentz transformation on the mass shell equation (15) in the comoving frame, rewritten 

as 2 4 2 2 

,,'2 ^2;/2 I / 2 i\, ,/2 ^ 

CO — CqK + [n — l)uj — "2 ■ [10) 

Note that the right hand side is a Lorentz scalar and the left hand side has a Lorentz 
invariant part cu''^ — c^k'^. The remaining part can be transformed to give 

- elk'' + (n^ - 1)7^ X (a; - k.vf = Hf^!^^ (74) 

where A;^ = —Icj is the wave four vector. There are several other ways in which 
Eq.(74) could be obtained, either by usual Legendrc transformations from Eq.(72) or by 
the Hamilton- Jacobi equation 9^''-§§:-§§7 = m^c^ with g^"' = $2 ^ [r^t^^ + (^2 - 1) V^'V]. 
We do not do it here. 

A further interesting result holds as a corollary to Sec. 4: For light waves in flat space, 
v'p = v'g = Cq = Vg, as expected. It should be noted that the Minkowski observers A and 
B can also be located in the asymptotic region and the entire analysis would remain the 
same. Prom the asymptotic vantage point, these observers can see that, near the horizon, 
n — > 00, then v'^, v'g —>■ 0, both for light and matter de Broglie waves. It is exactly here that 
we find that the conditions for optical black holes required by Leonhardt and Piwnicki 
[19] and Hau et al [20] are provided most naturally, that is, extremely low group velocity 
or high refractive index. In this respect, optical and gravitational black holes look indeed 
similar. Also, Vg = V , implying that, while A sees everything standstill at the horizon, B 
sees them moving away at the speed V because of his own relative motion. This is what 
we should really expect. 



cS-.^+(l-l).^(c.-^)^, (72, 



9 Summary and concluding remarks 

The present investigation is inspired by recent discoveries and analyses of light propa- 
gation in Bose- Einstein condensates [19,20]. The extremely low velocity of hght in such 
condensates lead to the possibihty of creating optical analogs of astrophysical black holes 
in the laboratory. In order to theoretically model this possibility, Leonhardt and Piwnicki 
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[6] proceed from the moving optical medium to an effective gravity field with a scalar 
curvature i? 7^ in which light propagation is shown to mimic that around a vortex core 
or optical black hole. Novello and SaUm [24] have recently shown that the propagation 
of photons in a nonlinear dielectric medium can also be described as a motion in an ef- 
fective spacetime geometry. Our approach here has been in the exact reverse direction: 
We proceed from the gravity field and arrive at an effective optical refractive medium 
and examine the theoretical consequences. The motion of this medium is caused by the 
relative motion between the observer B and the gravitating source. 

We must mention that works based on the above mentioned analogies provide some 
curious theoretical insights both in the real media and in the gravitational field, as a result 
of wisdom borrowed from one field and implanted into the other. This has been the basic 
philosophy of the present paper. Many more interesting results are known apart from the 
possibility of optical black holes stated above. For instance, an analysis in acoustic theory 
leads to a remarkable result that the Hawking radiation in black hole physics is not of 
dynamical, but kinematical origin (Visser, Ref.[9]). Conversely, a gravitational refractive 
index approach, similar in spirit to that of ours, has yielded the possibility of Cerenkov 
radiation in the outskirts of a wormhole throat [26-28]. In the present paper, we envisaged 
a nontrivial dispersive Fresnel drag coefficient in a gravity field. We must emphasize that 
these results are only of pedagogic interest at present. A further confirmation or otherwise 
of these results would establish the extent to which these analogies could actually be relied 
upon. 

We saw above how dispersion effects, both for massless and massive particles, appear 
naturally as a consequence of the systematic development of an effective medium approach 
to gravitational field. Various expressions for the drag coefficients result due to the use of 
VAL and different wavelengths used by the observer B. (See Refs.[l-3] for more detailed 
arguments on the question of the use of appropriate wavelength) . It is demonstrated that 
F2 is independent of A even in a dispersive medium for massive particles and that F2 
is identically zero. These results may have interesting implications for both optical and 
general relativity black holes. 

It does not seem easy to simulate real experiments, with our type of unbounded 
medium, that parallel those dealing with ordinary media like solid, fluid or superfluid. For 
this reason, we limited ourselves only to theoretical calculations of the drag coefficients 
and the expressions may be useful in the study of passage of light and cosmic particles in 
astrophysical media since what we actually see from the moving Earth is not what was 
originally sent from the source. This work is underway. 

We saw that the present analysis naturally complements the curved space analogy of a 
moving medium. Some of the key expressions in the comoving frame are indeed the same. 
Moreover, we can find a direct extension of the expressions to a genuine gravity field (Sec. 
8). The resulting Lagrangian and Hamiltonian describe the trajectories of a particle as 
viewed from the lab frame, say, a rocket. It also appears that the nomenclature "optical 
black hole" is quite apt as the conditions required for their creation are most naturally 
met near the gravitational horizon. This gives an indication that the behavior of the real 
optical medium should mimic that of our equivalent refractive medium around a coordi- 
nate singularity. A favorable situation is attained if light perceives the highly refractive 
real optical medium as dispersionless which, in our effective medium, is actually the case. 
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Leonhardt and Piwnicki [6] also make a similar statement in the context of their vortex 
analysis. It is interesting to note that an index of the form n — C/r, where C is a con- 
stant, when put in Eq.(4) yields orbits that resemble those around an optical vortex core 
[10]. A similar investigation with a different form of index has been reported also recently 
[29]. 
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